Phonon background versus analogue Hawking radiation in Bose-Einstein condensates 
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We determine the feasibility of detecting analogue Hawking radiation in a Bose-Einstein conden- 
sate in the presence of atom loss induced heating. We find that phonons created by three-body 
losses overshadow those due to analogue Hawking radiation. To overcome this problem, three-body 
losses may have to be suppressed, for example as proposed by Search et al. [Phys. Rev. Lett. 92 
140401 (2004)]. The reduction of losses to a few percent of their normal rate is typically sufficient to 
suppress the creation of loss phonons on the time scale of fast analogue Hawking phonon detection. 
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Introduction: An analogue of the Hawking effect [l|, 01 
should in principle exist in cold moving fluids [3|. A 
fluid whose velocity profile contains a transition from 
subsonic to supersonic flow will emit a thermal radi- 
ation of quantized sound waves (phonons). Gaseous 
Bose-Einstein condensates (BECs) were often consid- 
ered most promising for the observation of this phe- 
nomenon [4, 5, 6, 7, sj, [1, [T^. However, for accessible 
Hawking temperatures the condensate must have den- 
sities for which three-body losses become relevant 
These strongly constrain the time available for phonon 
detection already on the level of the mean-field. Here, 
we consider limitations arising from three-body loss in 
a quantum treatment (l^ and show that they are even 
more severe. 

The primary consequence of three-body loss is a re- 
duction of the condensate density in time, resulting in a 
decrease of the Hawking temperature. Three-body losses 
also heat the condensate by driving the many-body quan- 
tum state away from the Bogoliubov vacuum nj. This 
creates phonons, which contribute in general to a back- 
ground indistinguishable from analogue Hawking radia- 
tion Here we compare these two effects and show 
that in equilibrium the loss-phonons always overshadow 
those created by the Hawking effect. For the densities 
required to observe analogue Hawking radiation the 
loss-phonons are created on the same time-scale as the 
Hawking phonons. 

Our findings indicate that it may be necessary to sup- 
press three-body losses in a BEG in order to observe 
analogue Hawking radiation. Fortunately, suppression 
schemes exist [13, suppression will have a two- 

fold benefit: Firstly it increases the time-scale for the 
loss-heating to reach equilibrium, making it possible to 
conduct an experiment before they become relevant. Sec- 
ondly, we can employ higher density condensates and ob- 
tain a stronger Hawking signal, due to the reduced effect 
of loss on the mean-field. 

The surface at which the normal component of the con- 



densate flow exceeds the local speed of sound is termed 
sonic horizon As consequence of the causal discon- 
nection of the supersonic and subsonic regions, quantum 
field theory predicts particle creation [l7|. The ana- 
logue Hawking temperature that characterizes their ther- 
mal spectrum is given by Q 
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Here v is the flow speed, c the speed of sound, M = v/c 
the Mach number, Xh denotes a position on the horizon 
and n a normal vector to it. The underlying correspon- 
dence between the equations of motion for a scalar quan- 
tum fleld in curved space-time and for quantum phonons 
in a fluid can provide further analogies between cosmo- 
logical effects and phenomena in a fluid Q. 

Values of Th for typical fluids are very small. Hence 
Bose-Einstein condensates have been considered as prime 
candidates for an observation of the analogue Hawking 
effect, owing to their low temperatures. We have shown 
in Ref. 11] however, that to reach temperatures even 
of the order of 10 nK, condensates have to be typically 
driven into a regime where three-body recombination sig- 
nificantly affects the mean-field on a time scale of 50 ms. 
Here we employ Bogoliubov theory used in Ref. [11] to 
determine the consequences for the quantum field. 
Bogoliubov- de Gennes equations: We split the field oper- 
ator for bosonic atoms ^^(x) = '/'(x) -I- x(x) into a con- 
densate wavefunction (/)(x) — (^'(x)) and its quantum 
fluctuations x(x)- The condensate wave function obeys 
the Gross-Pitaevskii equation 

^Ccj^^ f-|^ + W + g\q^\A (3) 
ot \ 2m J 

where m is the atomic mass and g the interaction 
strength related to the scattering length as by g = 
Airh'^as/rn, while W denotes an external potential. We 
use / c?'^x](/)(x)]^ — A^cond, the number of condensate 
atoms. Where required, we use the notation p = jc^]^ 
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for the condensate density, v = ih/ (2m)[(V(/)*)(/>— </>* V0] 
for its velocity and c = y^gpjm for the speed of sound. 

We decompose the fluctuating component as x(x) = 
X]„[Mn(x)d„ -|-w*(x)aJ'J. The presence of a subscript dis- 
tinguishes references to the condensate velocity v from 
those to the mode v^. The functions Un,v„ 
the Bogoliubov-de Gennes (BdG) equations [l^, [l9|, |2C 
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BdG — 



c 



, (4) 



where Q = 1 — |(/))((/'|/A^cond projects onto the function 
space orthogonal to the condensate mode. The modes 
also must be normalized according to / d^x[|M„(x)p — 
|w„(x)p] = 1, for the fluctuations to obey bosonic com- 
mutation relations [d„,Q!j„] = 8n,m- 

For a homogeneous condensate in a quantization vol- 
ume V we have W = 4> = ji/g. The BdG equations 
are then solved by Uq{x) ~ -^u^e^^'^, Uq — ■^^=, 

Eq = yJW-<^c? + eg [Hi m. In the following we will 
assume that a sonic horizon is present in the condensate 
such that analogue Hawking radiation is created, but that 
the bulk condensate can still be considered as a homoge- 
neous reservoir. No details regarding how to achieve this 
situation are required here, but can be found in (llj . 
Loss induced phonons: Dziarmaga and Sacha have shown 
in Ref. [13] that besides a reduction of the condensate 
population, atom losses also result in creation of phonons 
since the many-body quantum state is driven away from 
the Bogoliubov vacuum |0), defined by q;„|0) = 0. These 
phonons will make a detection of analogue Hawking ra- 
diation more difficult. Thus the relative strength of the 
phonon sources has to be determined. We focus on three- 
body losses in what follows, since they are most promi- 
nent in BECs at high densities. 

Three-body recombination results in a molecule and 
an energetic atom. The excess energy due to molecu- 
lar binding, Ef, = h^/ma^, is split between the kinetic 
energies of molecule and fast atom in the ratio 1:2. 
When these kinetic energies suffice for both particles to 
leave the trap, their corresponding quantum fields can be 
eliminated from the picture [23[. One obtains an effec- 
tive master equation that describes the effect of the loss 
process on the quantum state of the remaining trapped 
atoms. The generalization for Z-body loss is [ij, [2^ 
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[H,p] d^V[¥]p, 
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where T>[a]p = dpa^ — d^dpjl — pa^djl. H is the usual 
Hamiltonian describing the conservative dynamics of the 
remaining trapped atoms. Importantly, 7; is the event- 
rate for a given loss process. Thus for example 73 = 
-f^s/S [l^l, where if 3 is the usual number loss rate for a 
condensate 1251. 



Inserting the expansion of ^ , Eq. ([5]) can be rewritten 
in terms of quasi-particle operators d^, ctln- C)ne then 
recognizes that each quasi-particle mode m is coupled 
to a heat reservoir and its occupation nm{t) will relax 
towards a thermal state For the high density con- 
densates that we consider, the three-body loss channel is 
strongly dominant (Tlj . If a single channel dominates, we 
can write Q3 



dn„r{t)/dt = -IjiaimN (t) [n,„(t) - ni 
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where N{t) is the number of condensate atoms. The 
coefficients aim and uim are determined by [l^ 
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with the condensate mode 0o defined hy 4) — N{t)(j)o. 
Eq. ([6]) evolves each occupation number towards the equi- 
librium value nim- This value itself is time dependent in 
general, but can be assumed to vary slowly if losses are 
not too strong [isj . 

In the following let the condensate be d-dimensional, 
with D = 3 — d tightly confined transverse dimensions. 
We still allow D = 0. The ideal setup for analogue 
Hawking radiation is not obvious: Elongated harmoni- 
cally trapped condensates in 3D complicate matters with 
a nontrivial transverse structure of the sonic horizon [Tlj , 
whereas quasi ID or 2D trapping must avoid a Tonks gas 
or quasi-condensate [1^ . This could conflict with tightly 
confining the required high densities, as shown later. 

In the longitudinal dimensions we imagine a homoge- 
nous condensate over a volume V. The three-dimensional 
vector X is decomposed as x — z + r±, where z is 
longitudinal and r_L transverse. With this splitting, 
we can write the condensate and Bogoliubov modes: 
00 = Ae-'-^L/^'^^/VV, Uq(x) = ?iqe*iMe-'-i/2'^VVl^, 
Wq(x) = Uqe*'i^Ae~'"i/^'^ /W. Here a is the ground 
state width of the transverse harmonic confinement and 
rj_ = Ir^l. We fix A by f d°r±A'^ exp[-rla^] = 1, 
which gives A ~ Ti-i^-d.'J/'^cr-i^-'i)/^ , After replacing the 
discrete index m by the continuous label q, Eq. ([7]) can 
be solved by 
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aig = 7r(3-'')(i-')/2l/i--f'?(''-3)/2cr(3-'')(i-'), (8) 
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Let us denote the "Z-body loss temperature" of this 
thermal state by Tu- We have nig = [e""/^""^" - 
hence ksTu = e,/ log -I- 1]. Since niq = \vq\'^ « 

mc/2hq « 1/^9 ^ 1 for phonons, we can finally ap- 
proximate ksTu = eqniq. Using Vq = a/ mc/{2hq) and 
eq = hqc for phononic wave numbers gives 
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The Hawking temperature is limited to [ll| 
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The factor S is defined by \WM\x=xh. ~ I'^Oi]^^^ where 
is the heahng length at the horizon [11]. We require 
S ^ 1 for hydrodynamic flow. We see that with both ef- 
fects in equilibrium the loss-temperature is always greater 
than Th- In a situation where atoms are continuously 
lost from the condensate, a rigorous equilibrium in which 
to interpret Eq. ^ does not exist. Nonetheless, when 
the loss is not too strong, we expect a quasi-equilibrium 
to apply. It is found that the actual heating can even 
shghtly exceed Eq. (H]) ll3|. 

Heating time scale: The equilibrium temperature asso- 
ciated with a loss process is independent of the corre- 
sponding loss (damping) rate, while the evolution to- 
wards equilibrium, Eq. is not [T3| . A harmonic os- 
cillator damped by a thermal bath behaves similar. In- 
serting the expression for aim one obtains 



dt 



= -l^^-'y'^fp[-Ut)M)~nim], (11) 



where we use the condensate density in reduced dimen- 
sions Pcond = and the effective loss rate 



atom 


"He 23Na «^Rb ^^^^Cs 


g X 10^" [Jm^] 
Ks X 10*2 [mVs] 

Pmax X 10-'^ 

Th [nK] 


15.7 1 0.5 0.66 
9000 2.12 32 130 
3.0 194 50 25 
3.8 16 2.2 1.4 


Vo [%] 


0.7 13 0.2 0.09 


hw±/kB [/iK] 
ajx/27r [kHz] 
Et/kB [iiK] 
Er/kB [fiK] 


0.34 1.4 0.2 0.25 
1.2 3 0.4 0.5 
25 125 16.6 10.4 

2100 2700 190 29 
9.5 22 2.1 0.68 



TABLE I; Comparison of common BEC species regarding 
analogue Hawking radiation and loss heating. The upper 
4 rows are reproduced from Ref. iT]. Nextly we indicate 
770, defined as fraction of the original three-body loss rate 
that allows Th = 30nK for At — 50 ms, while delaying 
the time scale for loss phonon production to r ~ 20 At, see 
Eq. (|17|l . The lower 5 rows show the hierarchy of energy 
scales mc^ <^ hLjj_ <^ mm{Ei,, Et) required to achieve weakly 
interacting quasi ID or 2D trapping, while still allowing three- 
body loss products to escape. The lo± axe selected to allow 
this. The mean field energies use p = pmax. 
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If we consider short time scales on which p^ondlO '^^^ 
treated as constant, the solution of Eq. (fTT|) is 



{t)^mm fl-exp(-?('^-i)/VPc;nV) 



(13) 



The time-scale on which the phonon population reaches 
its equilibrium value is therefore r = '^^~''-'^^7r''pJ(^d- 
Since the reduced-dimensional density is related to 
the three-dimensional peak density by 
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titles as 



(3-d)/2^-(3-d) 



P3D,poak — 

we can estimate r using 3D quan- 
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For usual condensate densities r can be quite large. One 
obtains r = 141s for ^^Na at psu ~ 10^" m'^ {d ^ 1). 
However, the densities required for reasonable analogue 
Hawking temperatures are significantly higher jlll] with 
according to Eq. pT|) much faster loss phonon creation. 
Let us parametrize the density as 



P - 



V'Sf/KsAt, 



(15) 



which implies that within a time At a fraction / of this 
density will be lost due to three-body recombination [ll| . 
We now choose At as the time which is required to de- 
tect the Hawking effect. We imply At ~ 50 ms in what 
follows, about the time required for quick spectroscopic 
phonon detection [13] . We arrive at a heating time scale 



under these conditions of 



T = At/3/. 



(16) 



This is of the order of the proposed measurement time 
At and hence too short, unless we have / ^ 1/3. How- 
ever as was found in [11], in the regime of such small 
loss the Hawking temperatures become problematically 
low. Hence, we usually have r ~ At. Alternatively we 
can calculate the time to create one phonon in a given 
spectral region, by the Hawking effect or by the loss. For 
the same parameter regime as above, one finds that the 
time scales are comparable. 

Suppression of three-body loss: It has been proposed to 
inhibit three-body loss processes in BECs by periodically 
hipping the phase of the weakly bound molecular state 
that causes the loss [l5| . The phase flip can be achieved 
using repeated 27r pulses of laser light resonant on an 
electronic exited state transition of the bound state. De- 
structive interference is then responsible for a reduction 
in three-body loss rates to only a few percent of their 
usual value. We now investigate whether this is suffi- 
cient to overcome three-body loss related obstacles to 
the creation of analogue Hawking radiation in BECs. 

By Eq. ([16]), we require a loss suppression that enables 
/ <C 1/3. We pick a specified target temperature and 
assume three-body loss was reduced to = rjK^, with 
< 77 < 1. The fractional loss within the measurement 
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time At will then be 

/ = [2Ti\kBTHf^^K^l^t/ig''] 77 = /or?. (17) 

For this we have eliminated p between Eq. (|10|1 and 
Eq. using c = gp/m. See table [T] for values of 
77 that allow analogue Hawking radiation at = 30 nK 
while separating the time-scales for phonon measurement 
and loss heating (r » At). 

A crucial indicator for the efficiency of the scheme pre- 
sented in [iHl is the number of laser pulses that fit into 
the average life-time of a molecule, before its quantum 
state is perturbed by a collision with a condensate atom. 
This life-time can be estimated as Tmoi = {i^p)~^, where 
K ~ lQ-^^w?s-^ for 87Rb [2I and ^SNa [29|. For densi- 
ties pRb = 5 X 10^°m~^, pNa = 2 X 10^-^771"^, we obtain 
Tmoi ^ lO^s. According to Ref. [31 this allows a loss 
reduction to a few percent. 

A side effect of the laser pulses that suppress three- 
body losses is one-body loss due to Rayleigh scattering 
of laser photons The increased one-body loss rate 

has been estimated as 71 ~ 0.1 s^"'^. We then see from 
Eq. that the time scale for one-body loss induced 
phonon creation is about 10 s and hence unproblematic. 
Reaction products of loss process: Given the importance 
of three-body losses, we have to address the evolution of 
the molecules and fast atoms created in the recombina- 
tion process. For Eq. ^ to be valid, it is required that 
they are energetic enough to leave the trap [23| . Also for 
the cosmological analogy to hold, we wish to avoid the 
complications of a coupled atom-molecular condensate. 



Finally, collisions between the loss products and remain- 
ing atoms would induce further unwanted heating if the 
loss products remained in the trap [30l |. 

To ensure that the molecules and fast atoms can leave 
the trap, we require the trap depth characterized by hw± 
to satisfy fiwi^ Ei,. To avoid trapping so tight that we 
enter the Tonks gas regime, the parameter 7 = [riidCid]"^ 
as to be much smaller than one [26]. Here = npa^ 
and — h/ y/mnidgid are the ID density and heal- 
ing length respectively, with a — {h/mLu±)~^/^ and 
<7id — .'7/(27rtT^). The condition can be reformulated as 
huj± ^ Et = fi^ y^Trp/2as/m. If we wish to study the 
sonic horizon in a quasi one or two dimensional setup, 
the strength of transverse confinement is constrained by 
mc^ <C huj± ^ min(i?f,, Et), where mc^ is the interaction 
energy of the confined condensate. Exemplary numbers 
for these energies are shown in tablelH which demonstrate 
that this hierarchy can usually only just be fulfilled. 
Conclusions: We have shown that loss induced phonons 
are an overwhelming background for analogue Hawking 
radiation in a Bose-Einstein condensate. To overcome 
this problem we suggest a moderate suppression of three- 
body losses. This can make the time-scale of loss induced 
phonon creation sufficiently long for a fast detection of 
the analogue Hawking effect. 
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